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Abstract 

Radiative decays of the 4> meson have been studied using a data sample of 
about 20 million 4> decays collected by the CMD-2 detector at VEPP-2M collider 
in Novosibirsk. From selected e + e~ — > ir + ir~j events the <p — * 7r + ir~j decay has 
been observed for the first time. 

Under the assumption that the intermediate fo (980)7 state dominates in the 
4> — > 7r + 7r~7 decay, the corresponding branching ratio is Br{4> — ► /o (980)7) = 
(1.93 ±0.46 ±0.50) x 10" 4 . 

Selected e + e~ — ► /i + ,u~7 events were used to obtain Br(cf> — ► /x + / u _ 7) = 
(1.43 ± 0.45 ± 0.14) x 10~ 5 for £ 7 > 20 MeV. 

Using the same data sample, upper limits at 90% CL have been obtained for 
the C-violating decay of the <fi: Br((f> — > pj) < 1.2 x 10~ 5 ; and for the P- and 
CP-violating decay of the ry: 

Br(r] -> tt+tt ) < 3.3 x 10~ 4 . 

Introduction 

The identification of scalar mesons and particularly the determination of the light- 
est scalar qq nonet is of extreme importance both for quark systematics and searches for 
non-gg objects like glueballs and multiquark states (see the discussion of this problem 
in the minireview on scalar mesons by S.Spanier and N.Tornqvist, p. 390 of Ref. [ ![). 
It is generally accepted that the probabilities of the electric dipole radiative transitions 
of the meson are crucial for the clarification of the nature of /o(980) and a (980) 
mesons [ @, 3 ■ A search for these decays has been earlier performed by the ND [ |J and 
CMD-2 [ H groups. A number of new results has been recently published by SND [ 
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§, and CMD-2 [ g | collaborations at the e+e" collider VEPP-2M [ [10[ where the 
radiative decays (f) — > 7r°7r°7, <fi — > ^7r 7 and — > 7r + 7r~7 have been observed for the 
first time. 

The mode with two charged pions has very large background because of the radiative 
processes e + e~ — > 7r + 7r~7 where a photon comes from initial electrons or from final 
pions. Therefore a signal from the /o (980)7 final state can be seen as an interference 

structure at the energy E 1 = — ^- — - ~ 40 MeV in the photon spectrum. The shape of 
this structure depends on the /o(980) mass and width. In our first paper [ [| a search for 
the 4> — > 7r + 7r~7 decay led to the upper limit of 3 x 10 -5 that looked at first puzzlingly 
low compared to the decay <p ~> 7r 7r°7 observed with the branching ratio of about 
1 x 10~ 4 [ |7J. As was shown in [ |], ||, [11], in case of the (p — ► / (980)7 intermediate 



state it could be explained by the destructive interference between bremsstrahlung 
processes and the decay. 

In this paper results of the study of the e + e~ — > 7r + 7r~7 process with the CMD-2 
detector [ |12| are presented. In total, the 14.2 pb~ l of data have been collected since 
1993 at 14 energy points around the <fi mass. For the analysis of the 7r + 7r~7 decay mode 
13.1 pb^ 1 corresponding to 19.7 x 10 6 (f> decays were used. Seven scans of this energy 
region were performed allowing control of systematic errors caused by possible detector 
instability. The results obtained from individual scans were found to be consistent. 



Selection of 7r + 7r 7 Events 

Event candidates were selected by requiring only two minimum ionizing tracks in 
the drift chamber (DC) and one or two photons with energy greater than 20 MeV in 
the Csl calorimeter. The following selection criteria were used: 

1. The average momentum of two charged particles is higher than 240 MeV/c to 
remove the background from K s — > 7r + 7r~ decays. 

2. Detected tracks have a polar angle between 1.05 and 2.1 radians so that they 
enter the inner muon system. 

3. The sum of the energy depositions of two clusters associated with two tracks is 
less than 450 MeV to remove Bhabha events. 

4. The radial distance of the closest approach of each track to the beam axis is less 
than 0.3 cm. 

5. The Z-coordinate of the vertex (along the beam) is within 10 cm from the detector 
center. This cut reduces cosmic ray background by a factor of two. 
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6. Detected photons have a polar angle between 0.85 and 2.25 radians so that they 
enter the "good" region in the Csl barrel calorimeter. This requirement sup- 
pressed the background from the photons emitted by initial electrons. 

7. Events with an invariant mass of two photons close to the vr° mass (|rn 77 — m^o | < 
40 MeV) were removed. 

After the above cuts the main background for the 7r + 7r _ 7 final state comes from: 
a) the radiative process e + e~ — > /i + /i _ 7, b) the decay <fi — > 7r + 7r~7r° when one of the 
photons from the 7r° escapes detection and c) collinear events e + e~ — > fi + fi~,7T + 7T~ in 
which secondary decays and interactions of muons or pions with the detector material 
produce a background cluster mimicking a photon. 
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Figure 1: Calorimeter response for events with one or no hits in the muon system (a) 
and for events selected as muons (b). 



Selection of 7 Events 

The inner muon system was used to separate muons from pions. The muon system 
uses streamer tubes grouped in two layers (inner and outer) with a 15 cm magnet yoke 
serving as an absorber and has 1-3 cm spatial resolution. The requirement of hits in 
the inner muon system for both charged particles selects muon events, together with 
some pion events in which both pions pass the calorimeter without nuclear interaction. 

Separation of pion and muon events in the Csl calorimeter is illustrated in Fig. [TJ, 
where scatter plots of the energy deposition of one track vs. that of the other one are 
presented for events with one or no hits in the muon system (a) and selected as muons 
(b). Energy depositions are corrected for the incident angle. Pions can have nuclear 
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interactions and in some cases leave more energy, while muons mostly exhibit dE/dx 
losses only. 
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Figure 2: a.Muon system efficiency vs. energy, b. Reconstruction efficiency for DC. 
c. Probability for two pions to be selected as muons. 



To select a cleaner sample of muon events, in addition to the information from 
the muon system both tracks were required to show only minimum ionizing energy 
deposition in the calorimeter (60-130 MeV). All the remaining events were considered as 
pion candidates. In this raw separation the pion sample N V7T1 contains muons because 
of some inefficiency of the muon system while the muon sample A^-y contains pions 
since a pion can reach the muon system without nuclear interaction. The observed 
numbers of events N n7ry and are related to the true numbers N® and A° M7 as: 
AjrTry = (1 — e M ) • iV° M7 + (1 — e w ) • N® nj , 

where e M and e n are the muon system efficiency and probability for a pion pair to be 
selected as muons respectively 

The magnitudes of e M and e n were determined by studying correlations between the 
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energy deposition in the Csl calorimeter and response of the muon system to collinear 
7r + 7r~ and fi + fi~ events. The results of this study are shown in Fig. |2] for one of the 
experimental scans. Figure 0a presents the muon system efficiency for different energy 
points. Only statistical errors are shown. The probability for two pions to be selected 
as muons is presented in Fig. |2]c and is independent of the possible instability of the 
detection system. 

The number of pions and muons has been corrected for the DC reconstruction 
efficiency shown in Fig. |]b. This efficiency was determined from collinear Bhabha 
events as described in [ |J . 

Our study of the efficiencies gives a correlated 3% systematic uncertainty in the 
final number of tt + tt~ and fi + [i~ events. 
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Figure 3: Simulated probability for pions to be selected as muons vs. photon energy. 



For events with photons e M and e n depend on photon energy. This dependence was 
studied using simulated 7r + 7r~7 events. The probability of two pions to be selected 
as muons varies with photon energy as presented in Fig. [| and was used to correct 
the final number of events with photons. At low photon energy this probability is 
consistent with that obtained from collinear events (see Fig. ^|c). 



Constrained Fit 

To reduce the background from collinear events as well as that from the three pion 
decay a constrained fit was used requiring total energy-momentum conservation for 
a three body decay. About 20% of the selected events had an additional photon. In 
this case the constrained fit was applied to both possible combinations and that with 
a minimum \ 2 was chosen. 
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An additional cut was applied to the photon direction: an azimuthal angle should 
be more than 0.25 radians from the charged track direction. This cut removed the 
remaining collinear events with a background photon which survived after the con- 
strained fit. 

The x 2 /d.f. distribution for events selected as muons had very small background 
and was found to be in good agreement with simulation as shown in Fig. |]a. Motivated 




X E 7 , MeV 

Figure 4: a. The x 2 /d.f. distribution for events selected as muons. The histogram is 
simulation, b. Photon spectrum for events selected as pions for x 2 /d.f.<3 (dark points) 
and for 3 < \ 2 /d-f- < 6 (open points). The histogram is simulation of background 
from 7T + 7r~7r events. 

by simulation, a cut that x 2 /d.f. be less than 3 was imposed for pion and muon events 
selecting 95% of signal events. 

After the above cut the pion sample still contained some background mainly from 
three pion decays. Such three pion background appears when one of the photons 
from the ir° has energy below 20 MeV and is not detected so that the event looks 
like a three body decay with the remaining photon energy higher than 150 MeV. The 
X 2 /d.f distribution for these events was flat and events with 3 < x 2 /d-f- < 6, shown 
in Fig. [|b by open points, were used to estimate the background spectrum. At each 
energy the background spectrum was subtracted from that for the signal candidates. 

From a constrained fit one can obtain an improved estimate for the photon energy. 
This effect was studied by using simulation and results are presented in Fig. |5]. Sim- 
ulation shows that after the constrained fit photons have energy resolution about 5 
MeV in the whole energy range instead of erg = 8% x E 1 Csl resolution as shown in 
Figs. [|a,b. Figure |5|c demonstrates the simulated photon detection efficiency in the 
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Csl calorimeter vs. photon energy. The overall detector efficiency vs. photon energy 
is shown in Fig. [5]d for 7r + 7r~7 events. 




Figure 5: Study of simulated 7r + 7r - 7 events, a. Difference between the photon energy 
measured in the calorimeter and the initially simulated one vs. photon energy, b. 
Difference between the kinematically reconstructed photon energy and the initially 
simulated one vs. photon energy, c. Simulated photon detection efficiency in the Csl 
calorimeter, d. Overall detector efficiency vs. photon energy. 



To extract the resonant contribution associated with the 0, two data sets were 
used. The first set ("0" region) for E c m from 1016.0 to 1023.2 MeV had the integrated 
luminosity of 9.24 pb^ 1 while the second one ("off-0" region) taken at E cm . = 996-1013 
and 1026-1060 MeV with the integrated luminosity of 3.89 pfe -1 containing less than 
3% (j) decays was used for a background estimate. 

Figure |6] presents photon spectra obtained after background subtraction and cor- 
rections for the detector efficiency at the "0" (a,c) and "off-0" (b,d) regions for pions 



and muons respectively. The solid line corresponds to the theoretical calculation [ ITT 
taking into account the integrated luminosity at each energy point and p — uo mixing 
for the bremsstrahlung process. A peak at 220 MeV is due to the radiation by initial 
electrons resulting in the process e + e~ — > p'-f, p — > 7r + 7r~. The photon energy range 
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Figure 6: a. Photon spectrum for tt + ti 7 events at the "0" region; Solid line is pure 



bremsstrahlung. Dashed lines include a possible <p signal according to [ 111, b. Photon 
spectrum for 7r + 7r~7 events at the "off-0" region; c,d. The same spectra for /i + /i - 7 
events. 



20 to 120 MeV has minimum background and some excess of events over the expected 
bremsstrahlung spectrum (the solid line in Fig. |6]a) can be seen in the 7r + 7r~7 sample 
at the "0" region. In total, 30175 7r + 7r~7 events and 27188 /i + /i~7 events have been 
selected in this energy range. 



Cross Section Study 

The cross section for each energy point was calculated as a = N ev /(L-e). N ev is 
the number of selected events with photons in the energy range from 20 to 120 MeV. 
The integrated luminosity L for each energy point was determined from Bhabha events 
with about 2% systematic accuracy [ ||. The detection efficiency e was obtained by 
simulation and the approximation shown by the solid line in Fig. |d was used to correct 
photon spectra. 

The resulting cross sections of the processes e + e~ — > it + -n~^ and e + e~ — > /i + /i~7 
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versus energy are presented in Fig. |7|a,b. Only statistical errors are shown. The 
systematic error in the experimental cross sections was estimated to be about 5% 
dominated by the uncertainty of the pion-muon separation efficiencies. 
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Figure 7: a. Cross section for e + e~ — > 7r + 7r~7. Lines are theoretical predictions in 
case of no direct decay (solid line) and best fit (dotted line), b. Cross section for 
e + e" — * /i + //~7 with the theoretical prediction. 

The cross sections can be described according to the calculations performed in j 



TTJI . It has been shown there that the amplitudes describing initial bremsstrahlung 



differ from those for the final bremsstrahlung and the 4> meson, so that the differential 
cross section for 7r + 7r~7 events is: 

« \A%(s } E y )\ 2 + \Al(s, E,) + M*, E v) ± e« ■ A fo (s, E 7 ) | 2 . 



Here A l £(s, E 7 ) and A(£(s, E y ) are the amplitudes describing initial and final bremsstrahlung 
processes, and the amplitude A$(s, E 7 ) introduces the influence of the upon the 
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photon propagator (a vacuum polarization term). This contribution gives rise to an 
interference pattern in the cross section at the mass and can be referred to as an 
"electromagnetic" decay 0^7^ 7r + 7r~7 or —>■ 7 — ► fi + fi~~f. The amplitude of the 
interference is determined by the meson leptonic coupling constant. 

The amplitude Af (s, E^) represents the possible decay into the 7r + 7r~7 final state 
via /07. For ^ + /x~7 events Af (s, E y ) = 0. The model considers the / meson as a two- 
or four-quark state or KK molecule depending on the values of the coupling constants 
g^/4vr and g^/47r. 

It also includes the interference of the amplitude with the final bremsstrahlung 
amplitude A[ r (s, E 7 ) as well as the correction for the 7r + 7r~ scattering in the final state [ 
|j~3| . The latter effect gives a contribution of about 15% to the branching ratio, but 
within our accuracy can be introduced as an additional phase shift ^ between the 
amplitudes of the and bremsstrahlung process. This shift was predicted to be 1.2-1.4 
radians, but the sign of the term was unknown and was determined from the fit. In 
general, all amplitudes have different dependence on s and E^. 

The initial bremsstrahlung process presented by the amplitude A % £(s, E y ) is sup- 
pressed by selecting photons transverse to the beam direction, but it still accounts for 
about 2/3 [ O of the observed e + e~ — > 7r + 7r~7 and one half of the e + e~ — > /i + /i - 7 
cross section. 

Integration over E 1 from 20 to 120 MeV gives a(s) which can be used for the cross 
section aproximation in Fig. |7|. 

Assuming no hadronic — > 7r + 7r~7 decay (Af (s, £' 7 )=0) a simple formula has been 
used where the final radiation amplitude of the bremsstrahlung process interferes with 
the Breit-Wigner amplitude from the vacuum polarization: 

*{*) = + • II - ^ ■ Af^\ 2 ; 
= s - m\ + iy/sT^s). 

Here s = 4E^ eamJ and are meson parameters, A ® is the decay amplitude 
in the peak and ip is the relative phase between the bremsstrahlung and decay 
amplitudes. 

This formula is valid under the natural assumption of the same photon spectrum 
for the final radiation amplitude of the bremsstrahlung process and the Breit-Wigner 
amplitude of the "electromagnetic" decay. 

The initial cr^f(s) and final u[ r (s) bremsstrahlung cross sections have different 
energy dependence for 20-120 MeV photons because of the pion formf actor energy 
behaviour. According to [ |TTJ the power functions cr£™(s) = 0.65 ■ cr ■ (wi^/^/s) 13 
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and cr( r (s) = 0.35 ■ cr • (m^/y/s) 9 can describe the bremsstrahlung process e + e~ — > 
7r + 7r _ 7 in our energy range. The parameter do represents the sum of initial and fi- 
nal bremsstrahlung cross section at the (f> mass. The function 0"£™(s) = &l r (s) = 
0.5 • 0"o(m^/s) was used for muons. 

The fit of experimental data with ctq, peak amplitude A$ and if) as free parameters 



shows good agreement of obtained cross sections with the theoretical calculations [ |TT 
al xp /af = 1.02 ± 0.02 ± 0.05, 
agv/of = 0.97 ± 0.02 ± 0.05 

for pions and muons respectively. The above values also indicate that the /i — n 
separation uncertainty does not exceed the estimated systematic error. 

Using the obtained amplitude one can calculate the peak cross section a ((f) — > 
7T7T7) = |A^°'| 2 • cr£,(m|) and the decay branching ratio can be calculated as Br((f> — > 
TTTf-f) = a ((f> — > 7fir^)/at ot , where af ot is the total peak cross section of the (f> resonance 
obtained from the leptonic width [ |J. The following results have been obtained after 
inserting corrections for angular acceptance equal to 0.68 and 0.42 for pions and muons 
respectively: 

Br((f) -> 7T+7T-7) = (0.41 ± 0.12 ± 0.04) x 10~ 4 , 

Br((p -> n+n~i) = (1.43 ± 0.45 ± 0.14) x 10~ 5 , 

if) = (0.46 ± 0.17) radians for pions and if) = (0.2 ± 0.4) radians for muons. 

The obtained values should be compared to the theoretical calculations based on 
"electromagnetic" (f> decays only [ [TT| : 

Br((p -> 7r+7T-7)=0.047xl0- 4 , 

Br((f) -> / u+/i"7)=l-15xl0- 5 . 
For muons the experimental value is in good agreement with theory, but for pions the 
value of the measured branching ratio is 9 times larger than the theoretical expectation 
and points to the presence of the hadronic decay of (f> to 7r + 7r~7. The obtained branch- 
ing ratios supersede our results in [ |[ based on the part of the total data sample. The 
corresponding fit curves are shown in Fig. [j]. 



Studies of Photon Spectra 

To search for the (f> — > /o (980)7 decay contribution to the observed <p — > 7r + 7r _ 7 
decay, photon energy spectra were studied. The photon spectra from the "0" region 
are shown in Fig. |3| for six c.m.energy points and photons in the 20-160 MeV energy 
range. The spectra were corrected for all experimental inefficiencies and normalized to 
the integrated luminosity. 

The signal from the decay of the (f> — > /o (980)7 is seen as a structure in the photon 
spectra at 40-60 MeV. Also shown by the dotted lines are the theoretical predictions 
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for the expected bremsstrahlung spectra [ [II, [13| including vacuum polarization . 

Because of the presence of the resonance in the photon spectra the branching ratio 
obtained in the previous chapter for the — > 7r + 7r~7 decay can be used only as some 
indication to the existence of the hadronic decay into this mode. 

It should be mentioned that the excess of events over the bremsstrahlung spectra 
cannot be used for the branching ratio calculation. As it was shown in [ ||, ||, 11 



the 



destructive interference with the bremsstrahlung process can reduce a visible signal in 
the selected energy range for photons. To demonstrate the effect of interference the 
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Figure 8: Photon spectra for six c.m.energy points in the "0" region. Lines are 
theoretical predictions for the four quark model (dashed line) and pure bremsstrahlung 
spectra (dotted line). 



photon spectra in Fig. || were fit as a group with the spectra calculated for each energy 
point using the differential cross section from [|2|, [Tl[ (see previous section). 



The number of events was not sufficient to keep all model parameters free and 
therefore the model parameters g^/47r and g^. 7r /47r were varied to keep the / (980) 
width at about 40 MeV [ [I[ . The fit had the following free parameters: the branching 
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ratio, /o(980) mass and phase shift The data can be described by the model only for 
the destructive interference between the fo amplitude and the bremsstrahlung process. 
The following results have been obtained: 

Br((f) -> / (980)7) = (1.93 ± 0.46 ± 0.50) x 10" 4 ; 

m fo = 976 ± 5 ± 6 MeV; 

* = 1.55 ±0.22 radians. 
The results above include a systematic error from the weak dependence on the /o(980) 
width. The relatively large value of the branching ratio obtained in the fit can not 
be accounted for by the models assuming the normal qq structure of / (980) [ [14| 
and is predicted only in case of its four quark structure. Note also that it significantly 
differs from the branching ratio obtained from the simple fit assuming the same photon 
spectra for all processes. 

Search for <f> — ► /cry Decay 

The selected 7r + 7r~7 events with photon energies from 100 to 300 Mev (see Fig. ^|) 
can be used to search for the C- violating decay 4> — > pj, p — > 7r + 7r~. The cross section 
vs. energy for the events with photons in the 100-300 MeV range is presented in Fig. 
|9|. The final state in the C-violating — > /ry decay has the same quantum numbers 
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Figure 9: <fi — > p7 search. The cross section for 7r + 7r _ 7 events with photons in the 
100-300 MeV range. The line is the best fit. 



as in the initial bremsstrahlung process which dominates in the cross section and to 
extract a possible signal the interference of these two processes is assumed: 
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<>%{*) 



ct(</>->P7) m^T^ p 
««(») ' I » 



where cr£"(s) = do • rn^/s. The comparison of the obtained cross section with the 
theoretical calculation [ [LTJ gives the ratio <To :rp /o"Q h =1.04±0.03. As a result of the fit, 
Br((j) -> p7) = (x(0 -> fn)loi ot = (0.3 ± 0.5) x 10~ 5 and ^ = -0.9 ± 1.0 radians have 
been obtained. Taking into account the 80% detection efficiency, the corresponding 
upper limit is: Br{<p -»• fry) < 1.2 X 10" 5 at 90% C.L. 

This result should be compared to the previous measurements which gave upper 
limits 7 x 10^ 4 [ § and 2 x 10~ 2 [ |XJ 



Search for 77 



7T + 7T" 



Decay 



The selected 7r + 7r~7 events can be used to search for the P- and CP- violating 
decay i] — > 7r + 7r~, where the rj comes from the radiative <fi — > 7/7 decay. From 19.7 
millions of decays used for the 7r + 7r - 7 channel analysis one could expect about 
248,000 events which decayed via the 777 channel. These P-, CP-violating decays 
should be observed as peaks in the invariant mass of two pions at m n7T = m^. Figure [10| 




1000 

eV 

Figure 10: Search for 77 — > 7r + 7r~ decay. The histogram in the box shows a simulated 
possible signal at 90% CL. 

represents experimental distributions of 7r + 7r~ masses from selected 7r + 7r~7 events. The 
line corresponds to a fit with a linear function and gaussian distribution representing 
a possible signal. It was found that the signal does not exceed 10 events at 90% CL. 
The histogram in the box shows a simulated signal from rj — ► 7r + 7r _ decays at 90% CL. 
The detection efficiency found by simulation was 0.124. The following result has been 
obtained: 
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Br(r/ -> 7T+7T-) < 3.3 x 10" 4 . 
which should be compared to the best previous limit of 9 x 1CT 4 [ 

Conclusions 

Using 13.1 pb~ 4 of data collected around the <p meson (about 20 millions of the <fi 
decays) e + e~ — > 7r + 7r~7 and e + e~ — ■> p + p~ r y events were selected. For the first time 
the decay — > 7r + 7r~7 has been observed in the 20-120 MeV photon energy range. The 
fit assuming that the <ft contributes only to the photon propagator (no direct <ft decay) 
gave the branching ratio: 

Br((j) -> 7T+7T-7) = (0.41 ± 0.12 ± 0.04) x 10~ 4 . 
This value is nine times higher than the expected and points to the presence of the 
hadronic decay of the (f) into this final state. The analysis of the photon spectra shows 
the presence of a resonance in the 7r + 7r~ system with a mass of about 980 MeV. The 
obtained branching ratio can be affected by the complicated interference of the hadronic 
4> decay with the bremsstrahlung process. The effect of interference was demonstrated 
by fitting photon spectra assuming the 4> — ► f (980)7 decay and the following results 
have been obtained: 

Br(<p -> / (980)7) = (1.93 ± 0.46 ± 0.50) x 10" 4 ; 

m fo = 976 ± 5 ± 6 MeV; 

* = 1.55 ± 0.22 radians, 
where \I/ is the relative phase of the amplitudes of the decay and the bremsstrahlung 
process taking into account the final state ir — ir interaction. 

According to the model [ 0, [11], 13], the obtained branching ratio can only be 



explained if / (980) is a four quark state. 
For muons the value 

Br((j) -> fi + fi--f) = (1.43 ± 0.45 ± 0.14) x 10~ 5 
has been obtained for the photon energies E 1 > 20MeV consistent with the theoretical 
expectations. 

For the C-violating decay of <fi — > p7 and P-, CP-violating decay of rj — ► 7T + ir~ the 
following upper limits at 90% CL have been obtained: 

Br((j) -> py) < 1.2 x 10" 5 , 

Br(j] -> 7T+7T-) < 3.3 x 10~ 4 . 
These results are the most stringent upper limits at the moment [ |I|]. 

Analysis of the / (980)7 intermediate state is continued in the following paper 
devoted to the studies of the and 777r°7 final states [ [TBI . 
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